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Drug delivery by micro- and nano-carriers enables controlled transport of pharmaceuticals to targeted sites. 
Even though carrier fabrication has made much progress recently, the delivery including controlled particle 
distribution and adhesion within the body remains a great challenge. The adhesion of carriers is strongly 
affected by their margination properties (migration toward walls) in the microvasculature. To investigate 
margination characteristics of carriers of different shapes and sizes and to elucidate the relevant physical 
mechanisms, we employ mesoscopic hydrodynamic simulations of blood flow. Particle margination is 
studied for a wide range of hematocrit values, vessel sizes, and flow rates, using two- and three-dimensional 
models. The simulations show that the margination properties of particles improve with increasing carrier 
size. Spherical particles yield slightly better margination than ellipsoidal carriers; however, ellipsoidal 
particles exhibit a slower rotational dynamics near a wall favoring their adhesion. In conclusion, 
micron-sized ellipsoidal particles are favorable for drug delivery in comparison with sub-micron spherical 
particles. 



The use of targeted micro- and nano-carriers for the delivery of imaging agents and drugs provides a 
promising strategy for early detection and treatment of diseases, e.g., of cancer^'^. However, the design of 
particles carrying different contrast agents and drugs as well as their physical delivery are very challenging 
tasks. Micro- and nano-particle fabrication, which needs to address several issues such as bio -compatibility, 
durability, binding to specific targets, and the ability of controlled release, has been strongly advanced in recent 
years^"^. Nevertheless, the development of efficient strategies for the delivery of carriers, including their distri- 
bution in the organism following systemic administration^ and their transport through biological barriers^" (e.g., 
microvascular walls, interstitial space, and cell membranes), requires a much more detailed understanding of the 
relevant physical and biological mechanisms^'^'^^'^^. 

Successful delivery of micro- and nano- carriers strongly depends on their efficient binding to specific targeted 
sites. Consequently, the distribution of carriers within vessel cross-sections plays an important role, since binding 
of carriers is only possible in case of direct particle-wall interactions. The cross-sectional distribution of micro- 
and nano-particles depends on several relevant parameters, which concern blood flow properties (such as flow 
rate, red blood cell deformability, and hematocrit - the volume fraction of red blood cells), vessel size, and particle 
characteristics (such as size, shape, and deformability). The migration of various suspended particles or cells 
toward walls in blood flow, which is often referred to as margination, has been observed experimentally for white 
blood cells^^'^^, platelets^^'^^, and rigid micro-particles^^'^^. Particle margination is mediated by red blood cells 
(RBCs), which migrate to the vessel center^^ due to hydrodynamic interactions with the walls (called lift force)^^'^^ 
leading to a RBC-free layer near the walls. More precisely, the occurrence of margination is a consequence of the 
competition between lift forces on RBCs and suspended particles, and their interactions in flow^^. However, the 
dependence of margination efficiency on particle size and shape remains largely unexplored so far. 

The role of particle size and shape in the efficient delivery is a multi-faceted problem. Large enough particles 
with a characteristic diameter (Dp) greater than about 4 jim may become trapped in the smallest capillaries of the 
body^^. In addition, recent experiments suggest that large particles with Dp > 3 |im are subject to an enhanced 
phagocytosis^^. However, recent microfluidic experiments^^ have shown that spheres with the size of 2 jam show a 
significantly higher adhesion density than particles with a size of 200 nm and 500 nm. Other experiments^^ 
indicate that liposomes with Dp < 70 nm and Dp > 300 nm have shorter circulation times than those having 
an intermediate size of Dp ~ 150-200 nm. Furthermore, nano-particles with a size below 20-30 nm are rapidly 
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excreted through the kidneys^^. Experiments with discoidal part- 
icles^^ have shown that they accumulate in the organs better than 
their spherical counterparts; however, particle internalization by 
macrophages appears to be worse for elongated particles^^. 
Adhesion of different particles has been studied experimentally^"'^^ 
and theoretically^^'^^ with the result that oblate ellipsoids are subject 
to stronger adhesion than spheres with the same volume. To better 
understand the adhesion potential of micro- and nano-particles, a 
quantitative description of particle margination under realistic blood 
flow conditions is required. 

In this work, we investigate the role of particle size and shape on 
the margination efficiency, and therefore on their adhesion potential. 
Several sizes ranging from about hundred nanometers to a few 
micrometers and two different shapes (spherical and ellipsoidal) 
are considered. The margination of micro- and nano-particles is 
studied numerically for a wide range of hematocrit values, vessel 
sizes, and flow rates using a combination of two-dimensional (2D) 
and three-dimensional (3D) models. Our results indicate that large 
particles possess a larger probability of being marginated than small 
particles. As the particle size becomes very small (less than about 
100-200 nm), the particle distribution within vessel cross-section 
can be described well by the plasma volume around flowing RBCs. 
Furthermore, spherical particles marginate better than ellipsoids, 
however the adhesion efficiency of ellipsoidal particles is expected 
to be superior in comparison to that of spheres due to their slower 
rotational dynamics near a wall. Finally, we discuss what size and 
shape of micro- and nano-carriers may be best suited for biomedical 
applications. 

Results 

Blood is modeled as a suspension of RBCs and micro- or nano- 
particles, while blood flow is studied in idealized microvessels using 
simulations in 2D and 3D, see Fig. 1 and Methods section for details. 
The 3D blood model has been shown to properly capture essential 
properties of blood flow in microchannels^^ as well as blood rheolo- 
gical characteristics^^'^^. To study micro- and nano-particle margina- 
tion for a wide range of conditions, we also exploit a 2D blood flow 
model due to its numerical efficiency; however, we will show that the 
2D model is able to qualitatively reproduce the required blood flow 
characteristics and the particle margination effect in comparison 
with the 3D model. 

Blood flow characteristics. The simulated system corresponds to a 
cylindrical micro vessel in 3D with the diameter W or to a channel in 
2D with the width W. We focus here on channel width W = 20 |im, 
but the cases ofW= 10 |im and W = 40 |im are also discussed. In 
flow direction, periodic boundary conditions are assumed and blood 
flow is driven by a constant force applied to all solvent particles, 
which is equivalent to a prescribed pressure drop. The hematocrit 




Figure 1 | Snapshots of cell and particle conformations in microchannels 

in 3D and 2D. RBCs are colored in red and suspended particles in blue, (a) 
3D simulation snapshot of blood flow for Ht = 0.3 and y*~39. (b) 2D 
simulation snapshot of blood flow for Hf = 0.3 and y*~29.3. 



Hf is defined as the volume fraction of RBCs. To characterize the flow 
strength, we define a non-dimensional shear rate in both 2D and 3D 
as 

y =yx = y , 1 

Kr 

where y = v/ W is the average shear rate (or pseudo shear rate) and v 
is the average flow velocity computed from the flow rate, while t 
defines a characteristic RBC relaxation time. Here, f\ is the 
solvent's dynamic viscosity, Dy = \/Ao/k in 3D and = Lq/tz in 
2D are the corresponding RBC diameters, where Aq is the RBC 
surface area in 3D and Lq is the cell contour length in 2D, and Kr is 
the RBC membrane bending rigidity. The RBCs are further 
characterized in 2D by the reduced area 

and in 3D by the reduced volume V* = 6Vo/ (ttDJ) =0.64, where 
Aq is enclosed RBC area in 2D and Vq is the enclosed RBC volume in 
3D. Typical values for healthy RBCs are = 6.5 jim in 3D, while 
= 6.1 |im in 2D, rj = 1.2 X 10"^ Pa s, and Kr lies within the range of 
SO-VOksT for the physiological temperature T = 37°C. Suspended 
micro- and nano-particles are characterized by the diameter for 
spheres and by the long axis for disks, denoted Dp in both cases. 

Particle margination in 2D and 3D. Margination of micro- and 
nano-particles in blood flow depends on H^, W, and y*. 
Figures 2(a),(b) illustrate the distribution of carriers of size Dp = 
O.lSDr (1.83 jam) for two Hf values in 3D. For better visibility, the 
carrier positions from a few snapshots are superimposed in the plot. 
The carrier surfaces are colored according to their radial position in 
the channel, with yellow color indicating a position near the channel 
center, while blue color corresponds to a position near the wall. 
Clearly, the carriers are marginating better for the case of larger Hf. 

Carrier positions in blood flow sampled over time lead to particle 
distributions, which reflect the probability of a particle to be at a 
certain distance from the wall. Figure 2(c) shows several center- of- 
mass distributions of circular particles in 2D with Dp = 0.3D^ 
(1.83 |im) for several Hf values and y*~29.3. The RBC-free layer 
(RBCFL) thickness, which is computed from simulation snapshots 
through the analysis of the RBC core boundary^'* similar to experi- 
mental measurements^^ (see Supplementary Fig. SI), is depicted by 
small arrows. The distributions have been averaged over the halves of 
the channel due to symmetry. Figure 2 shows that the carriers 
migrate into the RBCFL and remain quasi-trapped there. With 
increasing H^, the carriers marginate better, as indicated by the 
development of a strong peak in the distribution near the wall at 
y/W = 0, and the motion of the peak position towards the wall. 
This is due to a decrease in the RBCFL thickness leading to a smaller 
available space for the particles. This trend is in agreement with 
experimental observations^^ and simulations^^"^" of margination of 
blood platelets, which have a comparable size. 

To quantify and compare particle margination for a wide range of 
flow and particle parameters, we define the margination probability 
as a fraction of particles whose center- of-mass is located within the 
near- wall layer of thickness 3. The choice of 3 depends on the exact 
problem to be addressed, and several possibilities can be considered. 
To describe particle margination into the vicinity of a vessel wall, it is 
natural to select 3 to be the RBCFL thickness. Typical values of 
RBCFL thickness and their dependence on Hf are displayed in 
Supplementary Fig. S2. Figures 3(a), (b) present margination prob- 
ability diagrams of particles for a wide range of Hf and y* values 
corresponding to 3D and 2D simulations, respectively; the compar- 
ison shows that roughly 73^—1.2720. Particle margination strongly 
depends on as well as on shear rate. At low values, particle 
margination is expected to be weak, while at high Hf the margination 
might be also attenuated due to particle-RBC interactions near a wall. 
The latter effect has been described for a marginating white blood 
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Figure 2 | Particle distributions in blood flow. Illustrations of 3D simulations of blood flow for the shear rate of y* ~ 59 and different hematocrit values (a) 
Hf = 0.2 and (b) Ht = 0.4. RBCs are drawn in red, while spherical carriers with a size of Dp = O.ISD^ (1.83 [im) are colored according to their 
radial position r. For better contrast, carrier positions from several time instances are superimposed in the plots, (c) Center-of-mass distributions of 
carriers for various values at y*~29.3. 2D simulation results for circular particles with Dp = 0.3D^ (1.83 |am). The wall is at y/W = 0. The arrows 
indicate the boundary of the RBCFL for the different hematocrits, marked by corresponding colors. 



cefl^^ and is expected to subside for particles substantially smafler 
than a RBC, i.e. of sub -micrometer size. A pronounced dependence 
of particle margination on shear rate is observed at low flow rates. In 
the limit of very smafl flow rates (y* < 1), the RBC distribution 
should be nearly uniform, and therefore, the RBCFL and conse- 
quently particle margination should almost vanish. As the shear rate 
is increased, the RBCFL thickness grows rapidly^^, leading to a sub- 
stantial increase in particle margination. 

The simulated values of y* cover the range of flow rates character- 
istic for the venular part of microcirculation (y < 80 s~^ for W ~ 
20 jim), where it is estimated that y* < 90 in 3D {y* < 77 in 2D), 
while in arteriolar part the flow rates are higher {y* > 110 s~^ for 
W ~ 20 jim) with y* > 120 in 3D^^'^^. The considered range of shear 
rates is also relevant for tumor microvasculature, since blood flow 
velocities in tumors are much reduced in comparison to those under 
normal conditions, due to high geometric resistance and vessel per- 
meability^^'^^. Furthermore, the margination probability diagrams in 



Figs. 3(a),(b) show that the strongest particle margination occurs in 
the range of Hf = 0.25-0.6. This region has a considerable overlap 
with the characteristic hematocrits in the body's microvascular net- 
works in the range Hf = 0.2-0.4. A strong particle margination at 
high Hf values seems to be an advantage for drug delivery to tumors, 
since blood within tumor microvasculature is often subject to hemo- 
concentration due to plasma leakage^^. We also note that particle 
margination obtained from 3D simulations displays a higher mar- 
gination probability at lower values than that in the corresponding 
2D system. This difference arises from the variation in RBCFL thick- 
nesses in 2D and 3D systems. Thus, RBCFLs in 3D tubes are thinner 
than those in 2D channels for the same values due to cylindrical 
curvature of the geometry, which affects close-packing of flowing 
RBCs. In order to relate simulations with similar RBCFL thicknesses, 
3D margination data should be compared with 2D data at a larger 
hematocrit (by about 0.1-0.2), see Supplementary Fig. S2. In addi- 
tion, the 3D data also shows a decrease of particle margination at 




Figure 3 | Particle margination in 3D and 2D. Probability diagrams of particle margination with respect to y* and Ht in (a) 3D and (b) 2D, where the 
margination probability is defined as a probability of a particle center-of-mass to be within the RBCFL. The white squares ( □ ) indicate the values of Hf 
and y* for which simulations have been performed. 
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high shear rates. Nevertheless, both 2D and 3D simulations show 
qualitatively similar trends for the dependence of carrier margination 
on Ht and y* . Therefore, we conclude that 2D simulations are able to 
properly capture particle margination properties in blood flow. 
Further, we will often employ a 2D system due to its robustness 
and low computational cost in comparison with a 3D system. 

Dependence of margination on particle size. The discussion above 
considered the margination of micron- size particles. There is also a 
strong interest in nano-carriers, with sizes starting from several 
nanometers. Figures 4(a),(b) show margination diagrams of 
particles with = 0.15D^ (0.91 jim) and Dp = 0.04D^ (250 nm), 
respectively. The comparison of Figs. 4(a), (b) and Fig. 3(b) for Dp = 
0.3Dr (1.83 jim) reveals that the region of high margination 
probability becomes smaller with decreasing particle size. To 
illustrate the reason for the reduction in margination probability 
with decreasing particle size, we present in Fig. 4(c) the 
distributions of particles with different sizes for Hf = 0.3 and 
y*~29.3. For large enough particles, we observe a pronounced 
peak in the distribution next to the wall due to their interactions 
with RBCs, since their size is comparable with the RBCFL 
thickness. Even though small particles are also marginated, their 
distribution within the RBCFL is more uniform and their presence 
around the vessel center line is more probable than that for larger 
particles. Thus, the cumulative probability for a single particle to be 
within the RBCFL is lower for nano-carriers than that for micro- 
particles. Recent in vivo experiments^^ also support our numerical 
observations that particles with a size of about 1 jim are located 
closer to the vessel wall than smaller nano-particles. Noteworthy is 
that the distribution of the smallest particles with Dp = 0.04D^ closely 
approaches the distribution computed as the excess fluid volume of 
flowing RBCs. This indicates that the distribution of particles smaller 
in size than roughly 250 nm can be well approximated by the 



distribution of the blood plasma, and therefore, their margination 
properties can be directly inferred from local distributions. 

To decide on a suitable particle size for efficient drug delivery, a 
number of different considerations have to be taken into account. A 
direct interpretation of probabilities in Figs. 3(b) and 4(a), (b) sug- 
gests that larger particle sizes are more favorable for drug delivery 
due to their better margination properties. To further support this 
proposition, we consider another definition for the margination 
probability based on ^ = 0.5 Dp + s, which characterizes the fraction 
of carriers whose closest surface point is not further away from the 
wall than a distance s. We denote such a layer as "potential adhesion 
layer", since particle margination into a thin near-wall layer is a 
necessary precondition for adhesion. Even though the distance s is 
motivated by direct receptor-ligand interactions which occur within 
several nanometers, resolution restrictions in our mesoscale simu- 
lation approach do not allow the selection of smaller distances than 
approximately s = 0.03 ID^ in 3D or s = 0.033D^ in 2D, which 
corresponds to about 200 nm. Nevertheless, the distance of several 
hundred nanometers becomes relevant for particle-wall interactions 
in case of a carrier whose surface is decorated by tethered molecules^^. 
Another definition for margination probability can also be based on a 
fixed layer thickness S, thus it does not depend on or on particle 
size. Margination diagrams for this definition are illustrated in 
Supplementary Fig. S3. 

Figure 5(a) presents the margination probability into the potential 
adhesion layer (pj in 2D at y* — 29.3. At very small H^, the fraction of 
particles within the potential adhesion layer is small for all particle 
sizes; however, the smallest studied particles seem to be slightly more 
advantageous here. Remember that the interpretation of 2D mar- 
gination data for a RBCFL thickness with respect to the same RBCFL 
thickness in 3D requires a shift in values such that the range of 
Hf = 0.15-0.4 in 3D corresponds to approximately the range of Hf = 
0.3-0.6 in 2D. For the range of = 0.3-0.6, Fig. 5 clearly shows that 
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Figure 4 | Dependence of margination on particle size. Probability diagrams of particle margination in 2D for various Ht and y* values and for circular 
particles with the sizes (a) Dp = CISD^ (0.91 |am), (b) Dp = 0.04Dr (0.25 |am). The white squares ( □) indicate the values of Hf and y* for which 
simulation were performed. The margination probability is calculated based on the RBCFL thickness, (c) Distribution of particles with different sizes 
across the channel for Ht = 0.3 and y* —29.3. For small particles the distribution resembles the black solid curve computed as the blood-plasma volume. 
The arrow denotes position of the RBCFL boundary. 
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Figure 5 | Margination into a potential adhesion layer of thickness 200 nm. (a) Margination probability p^. The curves correspond to different particle 
sizes, where Dp = 0.63D^ (3.84 |am) is for an elliptic particle and the other curves are for circular particles. 2D simulation results for y*~29.3. (b-d) 
Margination diagrams for (b) Dp = 0.3D, (1.83 |am) (c) Dp = 0.15D, (0.91 |am), and (d) Dp = 0.04D, (0.25 |am). 



the fraction of large particles within the potential adhesion layer is 
much higher than that for small particles. The corresponding mar- 
gination diagrams are shown in Figs. 5 (b-d) and support the con- 
clusion that large particles marginate better for all considered shear 
rates. This indicates that micro -carriers are likely to be better for drug 
delivery than sub-micron particles. 

Dependence of margination on vessel size. To elucidate the effect of 
vessel diameter, we performed a number of simulations in 2D for two 
additional channel widths ( W = 10 jim and 40 jam) and two particle 
sizes {Dp = 0.1 5D^ and Dp = 0.3D^); corresponding simulation 
snapshots are displayed in the Supplementary Fig. S4. We consider 
margination into both the RBCFL and the potential adhesion layer. 

The pronounced dependence of particle margination properties 
on channel width for the potential adhesion layer is illustrated by a 
comparison of Fig. 5(b) and Fig. 6. For particles with a size of Dp = 
03Dr (1.83 |im), particle margination into the potential adhesion 
layer improves considerably as the channel size decreases due to 
the much smaller RBCFL thickness in narrow channels. Thus, par- 
ticle adhesion is expected to be more efficient in small vessels (i.e., 
capillaries) than in large vessels (i.e., venules and arterioles). 
Supplementary Fig. S5 supports this observation for particles with 
Dp = 0.1 5D^ (0.91 |im). Furthermore, a reduction of margination 
into the potential adhesion layer with decreasing particle size is 
found for all channel sizes. 



Particle margination based on the RBCFL thickness exhibits sim- 
ilar dependence on and flow rate for different channel widths, see 
Supplementary Fig. S6. For the channel width W = 10 jim, margina- 
tion into the RBCFL differs only slightly for different particle sizes up 
to Hf ~ 0.5. For Ht > 0.5, the particle radius might be larger than the 
RBCFL thickness, leading to an apparent decrease in margination for 
the large particles. For the cases W = 20 and 40 jim, where the 
RBCFL thickness is always larger than the particle radius, we observe 
that large particles marginate clearly better than small particles. 

Dependence of margination on particle shape. Advances in micro - 
and nano-particle fabrication facilitate the production of carriers of 
various shapes, including spherical, prolate and oblate ellipsoidal, 
and rod-like shapes^. However, advantages of different particle 
shapes for drug delivery are still to be explored. Thus, we 
investigate the effect of shape on the margination properties in 
blood flow. Figure 7 displays results of simulations in 2D for the 
margination probability (based on the RBCFL) of elliptic particles 
under various blood flow conditions in comparison to circular 
particles. The ellipse has an aspect ratio of about 7 and the longest 
diameter is Dp = 0.63D^ (3.84 |im); the enclosed area corresponds to 
the area of a circle with diameter Dp = 0.22D^ (1.35 |im). The plot 
indicates that margination of elliptic particles is slightly worse than 
that of circular particles. From these data we can also conclude that 
margination of the elliptic particles with a smaller aspect ratio than 7 
is similar to that presented in Fig. 7. However, since the largest 




Figure 6 | Margination for different channel widths. Margination into the potential adhesion layer based on ^ = 0.5Dp + 200 nm, for particles with size 
Dp = 0.3D, (1.83 |am) and two channel widths (a) W = 10 |am and (b) W = 40 \im. 
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Figure 7 | Dependence of margination on particle shape. Margination 
probabilities of ellipse-like particles (dashed lines) for various Hf and y* 
values in comparison to circular particles (solid lines) of the same area. The 
long axis of a 2D elliptic particle is Dp = 0.63Dr (3.84 |am) and the aspect 
ratio equals approximately 7. The margination probability is calculated 
based on the RBCFL thickness. 

diameter of the ellipse is larger than that of a circle with the same area, 
its margination into the potential adhesion layer, which is defined as 
a probability of a particle to be within a near-wall layer of thickness 
3 = 0.5Dp + 200 nm, appears to be considerably larger for ellipsoids 
than that for the corresponding sphere (see Fig. 5). We have also 
performed a number of 3D simulations with oblate ellipsoids having 
the same aspect ratio, which showed that their margination is 
qualitatively similar to that in 2D. 

Recent theoreticaP^'^^ and experimentaP°'^^ studies suggest that 
ellipsoidal particles possess better adhesion properties than spheres 
due to a larger contact area for adhesion interactions. It is also inter- 
esting to consider drag force on an ellipsoid or sphere in shear flow 
near a wall. In case of an ellipsoid close and parallel to a wall, the drag 
force is found to be smaller than that on a sphere with the same 
volume, which was estimated in separate simulations of a sphere 
and an ellipsoid in shear flow with fixed position. Thus, adhered 
ellipsoidal particles experience a lower drag force due to fluid flow 
than the corresponding spheres with the same volume. In conclu- 
sion, the current knowledge about adhesion of ellipsoidal particles 
and our simulation results on margination suggest that ellipsoidal 
particles are very likely a better choice for drug delivery than spher- 
ical particles. 

Dynamics of marginated particles. Local particle dynamics may 
also influence the margination and adhesion efficiency. Simula- 
tions in 3D show that the dynamics of marginated particles (i.e., 
within the RBCFL region) is different for spherical and ellipsoidal 
particles. A spherical particle is subject to a uniform rotation, while 
an ellipsoid displays tumbling dynamics. A quantitative analysis of 
the average angular velocities (w) of marginated particles shows that 
ellipsoidal particles rotate considerably slower within the RBCFL 



than spherical carriers, see Fig. 8. The comparison is made for a 
sphere and an ellipsoid of the same volume, while the long 
semiaxis of the ellipsoid is about twice the radius of the sphere. 
Within the RBCFL region it is plausible to assume a simple shear 
flow with the wall shear rate y^, which can be computed directly from 
a near-wall velocity profile or estimated from the pressure gradient 
applied to drive the flow. In Fig. 8, we compare the average angular 
velocities of marginated spherical and ellipsoidal particles with the 
theoretical predictions by Jeffery^^ for an oblate ellipsoidal particle in 
shear flow, which is given by 



where is the aspect ratio of major and minor axis. For a sphere, = 
1 which implies Ws = Y^/2, while for an ellipsoid, = 7 which results 
in CL)e~0.14y^, such that 0)^ < cOs- In addition, while the results for 
{oj) in Fig. 8(a) for a sphere are close to the theoretical results, the 
computed angular velocities for an ellipsoid in Fig. 8(b) are lower 
than the corresponding theoretical predictions due to the 
confinement of the ellipsoid between the wall and flowing RBCs. A 
lower rotational velocity of a particle leads to a longer interaction 
time between the particle and a wall. Thus, adhesion of ellipsoidal 
particles is expected to be more efficient than for spheres with a 
comparable size. In conclusion, a detailed analysis of dynamics of 
marginated particles further supports the proposition that ellipsoidal 
particles are likely to be better candidates for drug delivery. 

Discussion 

Particle margination in blood flow depends on particle size and 
shape, hematocrit, vessel size, and flow rate. Margination of spherical 
and ellipsoidal particles increases with increasing hematocrit, while 
their margination properties appear to be rather similar, where a 
sphere marginates slightly more efficient than an ellipsoid. The pre- 
sented diagrams show that larger particles have a higher margination 
probability in comparison to the smaller ones. Moreover, the distri- 
bution of very small particles with a diameter smaller than approxi- 
mately 250 nm is well represented by the blood plasma volume of 
RBCs. Margination of particles into the potential adhesion layer is 
found to be more pronounced in small vessels, indicating that par- 
ticle adhesion is likely to occur more often in capillaries than in 
arterioles and venules. 

The simulation results are in good qualitative agreement with 
several experimental observations^^'^^'^^'^°'^^'^^. For example, margina- 
tion of micro -particles has been observed to be more efficient than 
that of nano -particles in recent in vivo experiments^^. However, a 
detailed quantitative comparison is still difficult due to two reasons. 
On the one hand, the majority of the simulation results is obtained 
for 2D systems, which provide interesting insights into the relevant 
mechanisms, but have limited power for quantitative predictions for 
3D systems. On the other hand, experimental data on particle mar- 
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Figure 8 | Dynamics of marginated spherical and elliptical particles within the RBCFL. Comparison of average angular velocities (co) of (a) spherical and 
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gination in blood flow^^'^^ are very scarce and most of the available 
experimental investigations (e.g., Refs. 25, 30, 31) focus on carrier 
adhesion. Even though margination is a necessary pre-condition for 
particle adhesion to vessel walls, particle margination and adhesion 
are not equivalent, since carrier adhesion may also depend on other 
factors (e.g., specific targets, the receptor/ligand density and 
distribution). 

Clearly, the size and shape of drug carriers are important para- 
meters not only for margination, but also for their adhesion and 
further transport through biological barriers (e.g., internaUzation). 
Our simulations suggest that ellipsoidal particles are expected to 
adhere more efficiently than spherical carriers due to a larger surface 
for adhesive interactions and decelerated tumbling motion within 
the RBCFL. Therefore, our future numerical investigations will be 
focused on the adhesion ability of various particles in blood flow. 
Further requirements for efficient drug delivery include particle 
transport through vessel walls, interstitial space, and cell membranes. 
For instance, particle internalization by endothelial cells and intra- 
cellular trafficking have been shown to be most efficient for spherical 
sub-micron particles, rather than for micron-size carriers with an 
ellipsoidal shape^°. This observation points in the direction of smaller 
carrier to be most efficient for internalization. As a consequence, the 
concept of multi-stage drug-delivery carriers^'^, where a larger micro- 
particle incorporates a number of small nano-carriers, seems to be 
very promising. In this way, margination and carrier delivery or 
adhesion to a specific target within the microvasculature could be 
achieved using micro-particles, which would then be followed by the 
release of nano -particles into the tissue. In conclusion, tackling vari- 
ous drug- delivery challenges is a complex issue; its resolution 
requires an inter- disciplinary effort including in vitro and in vivo 
experiments and realistic numerical simulations. 

Methods 

Simulation method. We employ the dissipative particle dynamics (DPD) method^^'^^ 
for 2D simulations and the smoothed DPD (SDPD) method^^ for 3D simulations, 
where both methods are mesoscopic particle-based simulation approaches which 
properly capture hydrodynamics. Both simulated systems are represented by a 
collection of n point particles. The particles interact locally within a selected cutoff 
region through three pairwise forces denoted as conservative (F'^), dissipative (F^), 
and random (F^) forces. The time evolution of the velocity v, and position r, of particle 
/ with the mass m; is determined by the Newton's second law of motion drj = Vidt and 

dvi = — (Fp + Ff + Ff ) c/t. More information on the DPD formulation can be found 

in Refs. 51, 52, while the DPD parameters used in 2D simulations are presented in 
Supplementary Tab. 1. 

In the SDPD method, the forces are derived by a discretization of the Navier- Stokes 
equation similar to the smoothed particle hydrodynamics (SPH) method^"*, while the 
implementation of thermal fluctuations is analogous to that in DPD^^. The forces on 

particle i are given by the conservative force ¥f = + ojijTij, the 

dissipative force Ff = — yij(vij + rijrij-Vij), and the random force 

Ff = y^ ^. (Jij (^dWij + ^ [dM^ij] ^ 'iij. Thereby, cAVy is a matrix of independent 

Wiener increments and dWij is its traceless symmetric part. Here, pi andpj are particle 
pressures which are given by the equation of state p = po{pl poY — h withpo> Po> '^y and 
h being model parameters, see Supplementary Tab. S2. The particle density pi is 

calculated locally as p^ = W(^rij), where W{r) = Y^~~3" (^^ ~^ ^ ~~) (^^ ~ 

the Lucy function^"^ and is the cutoff radius. The weight function w(r) is determined 
by V W(r) = — rw(r). The coefficients and cTy define the strength of dissipative and 
random forces. The fluctuation-dissipation theorem to be satisfied requires 

Oij = 2JkBTyij. Finally, the friction coefficients are defined as = ^ '-^^ . 
V 3 PiPj 

Blood components. In 3D, a RBC membrane and suspended carriers are both 
modeled by a collection of discrete points, which are the vertices of a triangular 
network of springs on their membrane surface^^ The network assumes fixed 
connectivity with the potential energy defined as 

^membrane ^spring ~l~ ^bend ~l~ ^area ~l~ ^voU (3) 

which includes the spring's elastic energy Uspring> the bending energy Ui^endy and the 
area and volume conservation constraints Uarea and Uyoi- The spring forces mimic the 



elasticity of a membrane. The bending energy represents the bending resistance of a 
membrane, while the area and volume energies enforce area-incompressibility of a 
membrane and incompressibility of the inner cytosol, respectively. Detailed 
description of these potentials can be found in Ref 55, while all model parameters are 
given in Supplementary Tab. S3. 

In 2D, RBCs and micro- and nano-particles are modeled as closed bead-spring 
chains, which incorporate bending rigidity and an area constraint*\ The model 
parameters are presented in Supplementary Tab. S4. Carriers in 2D are modeled by a 
collection of particles (see Supplementary Tab. S5), which are constrained to 
maintain a rigid configuration. 

Simulation setup. The simulation setup consists of a single channel of cylindrical 
shape in 3D with diameter W = 20 |im and length of L = 12. 3D^. In 2D, a slit 
geometry with different widths W = 10, 20, and 40 |im and length L = 19. 5D^ 
(independent of W) is employed. The channel is filled with fluid particles and with N 
suspended carriers and ATrbc RBCs. The number of RBCs is computed according to 
channel hematocrit, which corresponds to the volume fraction of RBCs in 3D and to 
the area fraction of RBCs in 2D. The number of suspended particles for different 
simulations is provided in Supplementary Tab. S5. 

Boundary conditions. In the flow direction, periodic boundary conditions (BCs) 
were imposed, while in the other directions the suspension was confined by walls. The 
walls are modeled by frozen fluid particles with the same structure as the fluid, while 
the wall thickness is equal to r^. Thus, the interactions of fluid particles with wall 
particles are the same as the interactions between fluid particles, and the interactions 
of suspended carriers and cells with the wall are identical to those with a suspending 
fluid. To prevent wall penetration, fluid particles as well as vertices of RBCs and 
carriers are subject to reflection at the fluid- solid interface. We employed bounce- 
back reflections, because they provide a better approximation for the no-slip 
boundary conditions in comparison to specular reflection of particles. To ensure that 
no -slip boundary conditions are strictly satisfied, we also add a tangential adaptive 
shear force^^ which acts on the fluid particles in a near- wall layer of a thickness he = r^. 

Coupling between solvent and cells/carriers. Coupling between the fluid flow and 
cells/carriers is achieved through viscous friction^^ between cell vertices and the 
surrounding fluid particles, which is implemented via the DPD interactions and 
for both 2D and 3D simulations. The strength y of the dissipative force F^ for the 
interaction between a fluid particle and a membrane vertex is computed such that no- 
slip BCs are ensured. The derivation of y is based on the idealized case of linear shear 
flow over a flat part of a membrane with area A. In a continuum hydrodynamics 
description, the total shear force exerted by the fluid on the area A is equal to Arjy, 
where rj is the fluid's viscosity and y is the local wall shear-rate. The same fluid force 
has to be also transmitted onto a discrete membrane having vertices within the 
area A. The force on a single membrane vertex exerted by the sheared fluid can be 

found as Fv = ng{r)F^dV where n is the fluid number density, g{r) is the radial 

JVh 

distribution function of fluid particles with respect to the membrane particles, and 
is the half sphere volume of fluid above the membrane. Here, the total shear force on 
the area A is equal to N^Fy. The equality of AT^Fy = Arj^ results in an expression of the 
dissipative force coefficient in terms of the fluid density and viscosity, wall density 
Na/A, and r^. Under the assumption of linear shear flow the shear rate y cancels out. 
This formulation results in satisfaction of the no -slip BCs for the linear shear flow 
over a flat membrane; however, it also serves as an excellent approximation for no-slip 
at the membrane surface. Note that conservative interactions between fluid and 
membrane particles are turned off, which implies that the radial distribution function 
is structureless, g{r) = 1. In 2D, the surface area is replaced by a line of length L with 
Nl particles, and the half sphere volume is replaced by the half circle area Ay. 

Gathering statistics. The center-of-mass distributions of particles were calculated on 
the fly and were written to disk as sub-averages over short time intervals. Final 
averaging of the data is done during post-processing. To make sure that the final 
averaging of data starts from a time point which is sufficiently late for the system to be 
independent of the initial conditions, we have tested the sensitivity of final 
distributions to the choice of the starting time for averaging. 

Measuring RBC-free-layer thickness. To determine the RBC-free-layer (RBCFL) 
thickness, we measure the outer edge of the RBC core shown in Supplementary Fig. 
SI, which is similar to RBCFL measurements in experiments^^'^^. The data are 
averaged for many RBC snapshots at different times. In 3D, the RBC core edge is 
measured by projecting RBC vertices onto the x-y plane, where curves of the RBC core 
minimum and maximum are fitted similar to that in Supplementary Fig. SI. Here, we 
also perform averaging over different angular orientations (to exploit the cylindrical 
symmetry of the channel) in addition to the temporal averaging. 
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